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We consider the behavior of the tangential velocity of test particles moving in stable circular 
orbits in f(R) modified theories of gravity. A large number of observations at the galactic scale 
have shown that the rotational velocities of massive test particles (hydrogen clouds) tend towards 
constant values at large distances from the galactic center. We analyze the vacuum gravitational 
field equations in f(R) models in the constant velocity region, and the general form of the metric 
tensor is derived in a closed form. The resulting modification of the Einstein-Hilbert Lagrangian is 
of the form R +n , with the parameter n expressed in terms of the tangential velocity. Therefore we 
find that to explain the motion of test particles around galaxies requires only very mild deviations 
from classical general relativity, and that modified gravity can explain the galactic dynamics without 
the need of introducing dark matter. 
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I. INTRODUCTION 

Modern astrophysical and cosmological models are 
faced with two severe theoretical difficulties, that can 
be summarized as the dark energy and the dark mat- 
ter problems. Despite the fact that several suggestions 
have recently been proposed to overcome these issues, a 
satisfactory answer is yet to be obtained. However, in 
the context of dark matter, two observations, namely, 
the behavior of the galactic rotation curves and the mass 
discrepancy in galactic clusters, suggest the existence of a 
(non or weakly interacting) form of dark matter at galac- 
tic and extra-galactic scales. 

The galactic rotation curves of spiral galaxies [l| are 
probably the most striking evidences for the possible fail- 
ure of Newtonian gravity and of the general theory of 
relativity on galactic and intergalactic scales. In these 
galaxies, neutral hydrogen clouds are observed at large 
distances from the center, much beyond the extent of the 
luminous matter. As these clouds are moving in circular 
orbits with nearly constant tangential velocity i>t g , such 
orbits are maintained by the balance between the cen- 
trifugal acceleration v^ g /r and the gravitational attrac- 
tion GM(r)/r 2 of the total mass M{r) contained within 
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the radius r. This yields an expression for the galac- 
tic mass profile of the form M(r) = rv\JG, with the 
mass increasing linearly with r, even at large distances, 
where very little luminous matter has been detected 
This peculiar behavior of the rotation curves is usually 
explained by postulating the existence of dark matter, as- 
sumed to be a cold and pressureless medium, distributed 
in a spherical halo around the galaxies. 

There are many possible candidates for dark matter, 
the most popular ones being the weakly interacting mas- 
sive particles (WIMP). Their interaction cross sections 
with normal baryonic matter, although extremely small, 
are expected to be non-zero, and therefore it is believed 
to detect them directly [2j. However, no direct (non- 
gravitational) evidence for the existence of dark matter 
has been reported so far. It is important to emphasize 
that dark matter consisting of WIMP's may exist in the 
form of an Einstein cluster [3[ , or could possibly undergo 
a phase transition to form a Bose-Einstein condensate [J] . 
One cannot also a priori exclude the possibility that Ein- 
stein's (and Newton's) theory of gravity breaks down at 
galactic scales. In this context, several theoretical mod- 
els, based on a modification of Newton's law or of general 
relativity, have been proposed so far to explain the be- 
havior of the galactic rotation curves @, H, 0] ■ 

A promising avenue that has been extensively investi- 
gated recently are the f(R) modified theories of gravity, 
where the standard Einstein-Hilbert action is replaced 
by an arbitrary function of the Ricci scalar R Q ■ In this 
work we shall use the metric formalism, which consists 
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in varying the action with respect to g M „, although other 
alternative approaches have been considered in the liter- 
ature, namely, the Palatini formalism where the 
metric and the connections are treated as separate vari- 
ables; and the metric-affine formalism, where the mat- 
ter part of the action now depends and is varied with 
respect to the connection [loj . It has been suggested 
that these modified gravity models account for the late 
time acceleration of the universe [ll|, thus challenging 
the need for dark energy. However, the viability of the 
f{R) models proposed in the literature has been exten- 
sively analyzed [H, [H, [3 • I n this context, severe weak 
field constraints in the solar system ran ge seem to rule 
out most of the models proposed so far 1151 Hp , fl7l. fl8l| . 
although viable models do exist [H, [H, l2dl. l2l|. 

In order to be a viable theory, in addition to satisfying 
the solar system constraints, the proposed models should 
simultaneously account for the four distinct cosmolog- 
ical phases, namely, inflation, the radiation-dominated 
and matter-dominated epochs, and the late-time acceler- 
ated expansion |22j, and be consistent with cosmological 
structure formation observations [23], The issue of sta- 
bility [24[ also plays an important role in the viability 
of cosmological solutions [2(| HH, [Hj]. It is interesting 
to note that, recently, viable cosmological f(R) models 
were analyzed, and it was found that the latter models 
satisfying cosmological and local gravity constraints are 
practically indistinguishable from the ACDM model, at 
least at the background level [2l[ . 

The possibility that the galactic dynamics of massive 
test particles may be understood without the need for 
dark matter was also considered in the framework of f(R) 
modified theories of gravity [26|, [U [H[ . In the context 
of galactic dynamics, a version of f(R) gravity models 
admitting a modified Schwarzschild-de Sitter metric was 
analyzed in [2^ |. In the weak field limit one obtains a 
small logarithmic correction to the Newtonian potential, 
and a test star moving in such a spacetime acquires a 
constant asymptotic speed at large distances. It is inter- 
esting to note that the model has similar properties with 
MOND [5] . A model based on a generalized action with 
f(R) = R + R(R/Ro + 2/a)- 1 \n(R/R c ), where a, R 
and R c are constants, was proposed in [301 ]. In partic- 
ular, this model can describe the Pioneer anomaly and 
the flat rotation curves of the spiral galaxies. In a cos- 
mological context, the vacuum solution also results in a 
late time acceleration for the universe. The generaliza- 
tion of the virial theorem in f(R) modified gravity, using 
the collisionless Boltzmann equation, was considered in 

m. 

Within the framework of f(R) gravity, a model ex- 
hibiting an explicit coupling of an arbitrary function of 
R with the matter Lagrangian density was proposed re- 
cently 32]. Due to this coupling a connection between 
the problem of the rotation curve of galaxies, via a so- 
lution somewhat similar to the one put forward in the 
context of MOND, and the Pioneer anomaly is estab- 
lished. 



It is the purpose of the present paper to consider, from 
an exact analytic point of view, the problem of the galac- 
tic rotation curves in the framework of f(R) modified 
theories of gravity. In order to find an exact analytic 
description of the galactic dynamics of test particles in 
f(R) gravity models, we start from the general relativis- 
tic expression of the tangential velocity v tg of massive 
test particles in static and spherically symmetric space- 
times, moving in stable circular orbits around the galactic 
center. The rotational velocity is determined by the gu 
component of the metric tensor and of the radial distance 
only. 

We limit our analysis to the most important region 
for the galactic dynamics of test particles, namely, the 
region of constant rotational velocities. The constancy 
of the tangential velocity completely determines the form 
of gu, and consequently, the exact analytical form of the 
latter metric tensor component is completely determined 
from dynamical considerations. 

As a next step in our analysis of the geometry in the 
constant velocity region we consider the spherically sym- 
metric vacuum solutions of the gravitational field equa- 
tions in f(R) modified theories of gravity. By introducing 
several coordinate and functional transformations, the 
field equations can be reduced to an autonomous sys- 
tem of differential equations. By using the general form 
of g tt we obtain a second order differential equation fixing 
the functional form of g rr in the constant velocity region. 
This equation has as an exact solution g rr = constant, 
which gives g rr as a function of the observed tangential 
velocity only. The expressions of the Ricci scalar and of 
the function f(R) are also obtained. 

As a general conclusion of our study we find that to 
explain the flat galactic rotation curves, only small de- 
viations from standard general relativity are needed, so 
that f(R) oc R 1+V *s, a somewhat natural result in the 
context of modified gravity theories. 

As a possible observational test of our results we sug- 
gest the study of the lensing of light by galaxies in the 
constant velocity region. The study of lensing may, in 
principle, discriminate between the present model and 
other dark matter models. On the other hand, the deflec- 
tion angle in our model is of the same order of magnitude 
as the deflection angle in the standard isothermal sphere 
dark matter model, which is well tested observationally. 
This shows that the results obtained in this paper are 
consistent both theoretically and observationally. 

The present paper is organized as follows. The tan- 
gential velocity of a test particle in modified theories of 
gravity is derived in Section [TTJ The vacuum field equa- 
tions in the f(R) models are written down in Section [mi 
Two specific solutions of the field equations in the con- 
stant tangential velocity region are presented in Section 
IIV1 where some general properties of the basic equation 
describing the behavior of g rr are also discussed. We 
discuss and conclude our results in Section \V\ 
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II. THE MOTION OF MASSIVE TEST 
PARTICLES IN STABLE CIRCULAR ORBITS 

In order to obtain results which are relevant to the 
galactic dynamics, in the following, we restrict our study 
to the static and spherically symmetric metric given by 



ds 2 



^dt 2 + e x ^dr 2 +r 2 dn 2 7 



(1) 



where dfr 



de 2 



sin 8d(f) 2 . In the present paper, we 



use a system of units so that G = c = 1. The metric 
tensor coefficient e A< - r ) is constrained by the condition 
e A 0) > 1, Vr € [0, oo), while we assume that e v ^ satisfies 
the condition > 0, Vr € [0, oo). 

The Lagrangian C for a massive test particle reads 



C =2 



i {-e v i 2 + eV + r 2 f2 2 ) , 



(2) 



where the overdot denotes differentiation with respect 
to the affine parameter s. Since the metric tensor co- 
efficients do not explicitly depend on t and Q, the La- 
grangian ^ yields the following conserved quantities 
(generalized momenta) 33]: 



°r<- r H = E, 



2 n = L. 



(3) 



where E is related to the total energy of the particle and 
L to the total angular momentum. With the use of the 
conserved quantities, we obtain from Eq. ((2]) the geodesic 
equation for massive particles in the form 



e T r + e 



1 + K)=E 



(4) 



For the case of the motion of particles in circular and 
stable orbits the effective potential must satisfy the fol- 
lowing conditions: a) r = 0, representing circular mo- 
tion; b) dV e ff/dr = 0, providing extreme motion; c) 
d 2 V e ff/dr 2 \extr > 0, translating a stable orbit [33j |. 
Conditions a) and b) immediately provide the conserved 
quantities as 



Lr 

,-2 



and 



E 2 = e u 1 + -tt 



Lr 

,■2 



-T = -ire' V E 2 , 



(5) 



(6) 



respectively. Equivalently, these two equations can be 
rewritten as 



E 2 = 



1 - ru'/2 



L 2 = 



V/2 



1 - rz//2 



(7) 



We define the tangential velocity Vt g of a test particle, 
as measured in terms of the proper time [34| . that is, by 
an observer located at the given point, as 



"tg 



-u r 2 



dQ 
~dt 



L 2 

r?E 2 



(8) 



By using the constants of motion, we obtain the expres- 
sion of the tangential velocity of a test particle in a stable 
circular orbit [33| , given by 



TV 



(9) 



This simple expression which relates one of the metric 
components to the tangential velocity has three impor- 
tant properties. First of all, it is an exact general rel- 
ativistic expression valid for static and spherically sym- 
metric spacetimes. Secondly, it is interesting to note that 
the tangential velocity is sensitive to only one of the two 
metric functions, i.e., it is independent of the form of g rr . 
Lastly, since the motion of test particles is defined via the 
geodesic equations, this relation is independent of f(R). 
Even if we allow the modified theory of gravity to contain 
arbitrary contractions of the Ricci and Riemann tensors, 
the above equations would still hold exactly. 

In regions with constant tangential velocity the metric 
function v is fixed by the condition vt g ~ constant, and 
by integrating Eq. ([9]), we find 



v = 2v 2 g In 



(10) 



where ro is a constant of integration. Therefore the most 
general static and spherically symmetric metric in the 
constant tangential velocity regions can be written as 



ds 2 



dt 2 + e Mr) dr 2 +r 2 dfl 2 . 



(11) 



Note, however, that the galactic rotation curves gen- 
erally show a more complicated dynamics 1]. Therefore, 
to obtain a more accurate and realistic description of 
the particle's motion, which can be extended beyond the 
constant velocity region, more general expressions of the 
metric coefficients are needed. However, in the present 
paper we restrict our analysis to the constant velocity 
region, which provides interesting results. 

In order to test the consistency of our results, we con- 
sider the Newtonian limit of the model. The assumption 
of small velocities of the particles requires that the grav- 
itational field be weak. In the Newtonian limit the gu 
component of the metric tensor is given by e v « l + 2<i>Ar, 
where is the Newtonian gravitational potential sat- 
isfying the Poisson equation A<f>Ar = 4irp [34[- In the 
constant velocity region the mass M(r) of the dark mat- 
ter and the energy density p vary with the distance as 
M{r) = v 2 g r and p — v 2 g /Aitr 2 , respectively. Therefore 
the Poisson equation is given by 



J_d_ 

r 2 dr 



,d<P 



dr 



u _tg 



and has the general solution 
$jv(r) 



2 , r C N 
u tg ln • 

6 tq r 



(12) 



(13) 
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where Cn and tq are arbitrary constants of integration. 
In the limit of large r, corresponding to the constant 
velocity regions around galaxies, the Newtonian potential 
is given by 



/ = Fe 



v" -\{v' + X')v' --X' +(1/ + -\ 
2 r \ r J 



~F 



(20) 



$jv(r) 



V tg In 



ro 



(14) 



reflecting a logarithmic dependence on the radial distance 
r. On the other hand, the gu component of the metric 
tensor can be represented in the constant velocity "dark 
matter" region as 



e- w — 



1 



2< ln 



exp 

r 
ro 



In [ - 

ro 



= l + 2<S> N (r). 



(15) 



Therefore the model has a well-defined Newtonian 
limit, and the metric given by Eq. (|11|) can indeed be 
used to describe the geometry of the spacetime in the 
dark matter dominated regions. 



III. VACUUM FIELD EQUATIONS IN f(R) 
GRAVITY 

The action for the modified theories of gravity consid- 
ered in this work takes the following form 



S 



f(R)V=9 d 4 x, 



(16) 



where f(R) is an arbitrary analytical function of the Ricci 
scalar R. Note that we are only interested in the vac- 
uum case, and therefore we have not added a matter 
Lagrangian to the action. 

Varying the action with respect to the metric g^ yields 
the following field equations 

F(R)R^ - Lf(R)g^ - (V^V, - g^D) F(R) = 0, 

(17) 

where we have denoted F(R) — df(R)/dR. Note that 
the covariant derivative of these field equations vanishes 
for all f(R) by means of the generalized Bianchi identi- 
ties [32l. [351] . 

For a static and spherically symmetric metric of the 
form given by Eq. |T]), the field equations of the f(R) 
gravity in vacuum can be expressed as (29l . [36j 



F"--{u' + \')F'- {v ' + XI) F = Q, 
2 r 



(18) 



12 



-^' + A') 



= -2- 



2 

r 

El 



+ A' + 



EL 

~F : 



(19) 



F" 

I! = 2- 3e- A \ — 



It is useful to introduce a new variable 77 by means of 
the following transformation 



?/ = In r. 



(22) 



Therefore, the field equations Eqs. (fT8 |) -([2T |) take the 
form 



d 2 F 
dr/ 2 



1 + r 



dv 


dXY 


dF 


(dv 


dX 


dr\ 


dr)J 


df] 


\dri 




dr/ 



F = 0, 
(23) 



d 2 v dv 1 / dv dX\(dv 

dr] 2 dr] \drj J 2 \ drj dr] J \dr] 



+ 2 1 



x )=-2 



f = Fe- X - 271 



fv 
dr) 2 



1 d 2 F 
F~drf 



dv 1 

drj ~ 2 



dX + \ 1 dF 
dr] J F dr] ' 



dv dX\ dv 
dr] dr/ J drj 



R 



3e 



dr] 



1 d 2 F 



dv 
dq 



1 dF 
~F~drJ 



(25) 



F dr] 2 

1 (dv dX 

2 \dr] dr] 



1 dF 
F~drJ 



1 dF 
F~drJ 



(26) 



As a result of introducing the new variable the basic 
field equations Eqs. (|23p and (j2~4")) are independent of the 
radial coordinate ry. It is also very useful to introduce a 
formal representation of the function F as 



F (v) = F o exp 



u{r])dri 



(27) 



where u is a new function of 77, and Fo is an arbitrary 
constant, so that (l/F)dF/dr) = u, (l/F)d 2 F/drj 2 = 
du/dr] + u 2 . Hence Eq. (|23|) can be written as 



du 
dq 



u 2 ~ 



dv dX\ 
dr] drj J 



dv dX 
drj drj 



= 0. 
(28) 



This equation is a Riccati type first order differential 
equation. To find its general solution the knowledge of a 
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particular solution is required. By using the function u, 
Eq. can be written as 



dv ( dv 



1 / dv d\\ I dv 



drj 2 drj ~^ \drj J 2 \drj ^~ dr] ) \dr] 



du 



+ 2 (1 - e A ) = -2— - 2u 2 + — + 4 

" ' rirt-t \ /-fry-, 



drj 



dX 



dij 



(29) 



Substituting the term du/djj + u 2 in Eq. (|29]) . with the 
use of Eq. (|2"5|) , we obtain 



1 ^ dv\ , A , d 2 i^ / dv 

2 drj J ' d?7 2 drj \ drj 



dv dX^ 
drj drj , 



1 (iz/\ 



The general solution of Eqs. (|28|) and (f30|) provides 
the general static and vacuum solution of the modified 
field equations for the case of the f(R) gravity models. 
Once v{rj) and X(rj) are specified, one can immediately 
obtain u and consequently (by integration) F, as well as 
all the other relevant physical quantities. If the function 
F and the metric tensor coefficients are known, / can be 
obtained as a function of R from Eqs. (f2"5| and (|2l)|) in a 
parametric form, as f — f{i]), R — R(jj). 



IV. DARK MATTER AS A GEOMETRIC 
EFFECT IN f(R) MODELS 



equation 



1 



ld?X 

2 drj 2 1 - m \ 2 



3 \ 2 \ dX 

e — to + 1 — m ] — 
drj 



1 



(1-mY 



1-e' 



, 2 3m 2 - 1 



1-e" 



(1-to) z 
m 2 (2m 2 — l) 

— ^ >-2m = 0, (35) 

(1 — m) 



which must be satisfied by the metric coefficient A (77) 
in the "dark matter" dominated regions in f(R) gravity 
models. 

From astrophysical observations it is known that the 
tangential velocity of test particles in circular stable or- 
bits around the galactic center is of the order of i>t g ~ 
200 - 300 km/s [1]. Hence, we have m = v\ « 10~ 6 , 
and all the terms containing m 2 and to 4 in Eq. (|35|) can 
be neglected within a very good approximation. Consid- 
ering, thus, a first order approximation in m, Eq. (|35|) 
describing the metric coefficient exp(A) reduces to 



ld 2 A 1 



2 dry 2 1 - m \ 2 
x(l-e*) S 



, dX 

-e A + 1 -m — 



(1 - m)' 



dr/ (l - m ) 



(1 - e A ) - 2m = 0. (36) 



In the following analysis we will consider two classes 
of solutions of Eq. (|36|) . 



In order to obtain a geometric interpretation of dark 
matter in f(R) gravity, we start with the metric coeffi- 
cients v and A in the dark matter dominated region at 
the galactic scale, given by Eq. (fTTj) . which, using the 
variable 77, can be represented in the following form 

v = 2m (lnr/ — In 770) , A = A (77) , (31) 

where 770 = hiro = constant, and for notational simplic- 
ity we have denoted 



m = v t = constant. 



(32) 



For this form of the metric Eqs. (|28[) and ([3H]) become 



du 



drj 



1 dX 



dX 



l + m + --- ) u-2m- ^ = 0, (33) 
2 drj J drj 



and 



2 (1 - to) u = 2 (1 - e A ) + 2m 2 + (1 - to) — , (34) 

drj 

respectively. By taking the derivative with respect to 77 
of Eq. (|34p and substituting the resulting du/ dij and u in 
Eq. (f3"3"]l provides the following second order differential 



A. The case of the constant X 

Equation (|36[) admits an exact solution of the form 
A = constant. By denoting 1 — e A = 5, it follows that 
considering a first order approximation in to, then 5 sat- 
isfies the following second order algebraic equation 



6 - 2m = 0, 



(37) 



which has only one physical solution, namely, <5 = — 2to 
(the other solution contradicts the condition e A > 1, 
which essentially represents the positivity of the mass). 
Therefore, in the first order approximation the metric 
coefficient e A in the "dark matter" region becomes 



1 



2<, 



1 



2<. 



(38) 



Equation (|38p has a straightforward physical interpre- 
tation. Since in the Newtonian approximation v 2 g = 
GM(r)/r, where M(r) is the total mass of the galaxy, 
we obtain the metric in a form which is very simi- 
lar to the Schwarzschild solution of general relativity, 
i.e., e~ A w 1 — 2GM(r)/r. On the other hand, since 
v 2 g = constant, the mass within the radius r must in- 
crease so that M(r) ~ r. This "mass", which is linearly 
increasing with the distance, is usually interpreted as due 
to the presence of the dark matter. 



6 



In the present approach, the effect of the appearance 
of a dark "mass" is of a purely geometric origin, result- 
ing from the modification of the basic equations of the 
gravitational field. Nevertheless, one can formally define 
a mass M(r) for the dark matter, despite the fact that 
the origin of this "mass" cannot be related to physical 
particles. Thus, the rotational galactic curves can be 
naturally explained in f(R) gravity models without in- 
troducing any additional hypothesis. The galaxy is em- 
bedded in a modified spherically symmetric geometry, 
generated by the non-zero contributions of the modified 
gravitational action. The extra-terms act as a "matter" 
distribution outside the galaxy. 

Once the metric coefficients are known, Eq. (|33p pro- 
vides the function 14(77) to first order in m as 



u{rj) 



1 + m vT+IOto 



tanh 



\J\ + 10m (77 — 771 ) 



(39) 



where 771 is an arbitrary integration constant. The func- 
tion -F(t7) can be found as 



F (v) = F o exp 



1 



77 I cosh 



VI + 10m (77 - 771) 



(40) 



For large r, by taking one arbitrary integration con- 
stant equal to zero, without a significant loss of gener- 
ality, so that we keep only the decreasing term in the 
expression of F, we obtain 



(41) 



F{r) = C x r~ 



where C\ is an arbitrary constant. Finally, from Eqs. l[20l 
and (1211) we find 



f(R) = f R 1+m , 
where f is a constant. 



(42) 



B. First order corrections and general properties 

of A 



In order to obtain a better description of the astro- 
physical observations at the galactic scale, and also tak- 
ing into account the discussion of the previous Section, 
we may also consider a first order correction of the met- 
ric coefficient exp(A), by assuming that A is small. Hence 
we may approximate the term exp(A) as exp(A) ~ 1 + A. 
where A< 1. Moreover, the condition 5/2^>3A/2 — m 
is also satisfied. Therefore, by neglecting the quadratic 
term A~, Eq. (|36[) takes the form of an ordinary linear 
second order inhomogeneous differential equation, which 
can be written as 

(1 _ m f fL£ _ 5(1 _ m )^ + 2 A- 4777 = 0, (43) 



dt] 



and with the general solution given by 

A (77) = 2777 + C+ exp (s+77) + C— exp [s—T]) , (44) 
where C± are arbitrary constants of integration, and 

5±VT7 



2(1 - m)' 



(45) 



Hence, to first order in both m and A, the metric tensor 
coefficient exp (A) can be represented as 



e A = 1 + 2m + C*+r s + + C^r s ~ . 



(46) 



However, in this case the Riccati equation for u, 
Eq. l[2"g]). cannot be solved exactly, and numerical meth- 
ods are needed to further investigate the physical prop- 
erties of this solution. 

By introducing a new variable v = dX/drj, then 
Eq. (|36p can be transformed to a first order differential 
equation of the form 



1 dv 1 



2 dX 1 - m V 2 



— e + 1 — m ] v 



1 



(1 — my 



x (1-eT 



(1 - m)' 



(1 - e A ) - 2m = 0. (47) 



With the help of the transformation w = 1/v and by 
introducing a new independent variable — exp (A), we 
obtain 



dw 
~dB 



1 — 777 



1 



1 



(1-ro)' 



x(l 



(1 - my 



(1-9)- 2m 



w 3 = 0. 



(48) 



Hence, Eq. (136[) has been transformed to a first order 
nonlinear second kind Abel differential equation of the 
form dw/dO = A(9)w 3 +B(9)w 2 [37j]. As it is known from 
the theory of the Abel differential equations, an equation 
of this form has an exact solution if and only if the con- 
dition d [A(6UB(d)\ /d6 = kB{6) is satisfied, where k is 
a constant [37] ■ A simple calculation using the explicit 
forms of the functions A(0) and B(9) corresponding to 
Eq. (|48[) shows that this condition cannot be satisfied for 
all values of m. 

Therefore, it follows that the only exact solution of 
Eq. (|36[) describing the behavior of the metric coefficient 
exp(A) in the dark matter region in f[R) gravity models 
to first order in the tangential velocity is A = constant. 
To fully investigate the general behavior of this equation 
one must use numerical methods. 



V. DISCUSSIONS AND FINAL REMARKS 

From astrophysical observations and from their in- 
terpretation in the framework of the phenomenological 
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Modified Newtonian Dynamics (MOND) approach p|, 
it is known that the acceleration needed to explain the 
observed rotation curves is of the order of 10~ 10 m/s 2 , 
which can be regarded as only a small deviation from 
general relativity. Likewise, the Pioneer anomaly, which 
is also of the order of 10 -10 m/s 2 , would imply only small 
modifications of gravity. 

On the other hand, recent laboratory tests have con- 
firmed that Newton's second law is in "good agreement" 
with accelerations of the order of 5 x 1CP 14 m/s 2 [38|. 
Similar constraints have also been obtained for the in- 
verse square law, where it was shown that Newton's law 
holds down to a length scale of 56 ^m 39]. Hence, it fol- 
lows from these observations and experiments that in the 
R n modified theories of gravity the parameter n should 
be of the form 1 + e with eCl. 

In the present paper, we have considered the gravita- 
tional field equations in f(R) modified theories of gravi- 
tation in the flat rotation curves region. Although this is 
a rather strong assumption, it is valid for at least a sig- 
nificant region of the total velocity profile. As expected, 
the deviations from standard general relativity are only 
mild and the exponent of the power law modified gravity 
takes the form n = 1 + v 2 g . With the use of Eq. (|42f we 
obtain that this theory is defined by the action 

S= I ' f R l+ <^d l x, (49) 

where /o is a positive constant, given in terms of the tan- 
gential velocity, which can be obtained by using Eqs. (|20|) 
and (f2"Tjl . and v 2 g is a number of the order of 10 -6 . The 
function f(R) for this modified gravity model can also be 
approximated as f(R) = f R l+Vi * » /o-R(l + u 2 g lni?). 
Hence the correction terms to the standard Einstein- 
Hilbert action is logarithmic. 

The weak field limit of the f(R) generalized gravity 
models has been discussed recently, for star-like objects, 
in [ItJ and [Hj], respectively. By assuming that f(R) 
is an analytical function at the constant curvature i?o, 
that TTitpv -C 1, where is the effective mass of the 
scalar degree of freedom of the theory, and that the 
fluid is pressureless, the post-Newtonian potentials ^ (r) 
and $ (r) are obtained for a metric of the form ds 2 = 
- [1 - 2* (r)] dt 2 + [1 + 2$(r)] dr 2 +r 2 (d0 2 + sin 2 Ody 2 ) . 
One may then find the behavior of W (r) and $(r) out- 
side the star. The analysis leads to a value of 7 = 1/2 
for the Post-Newtonian parameter 7, which from So- 
lar System observations is known to have a value of 
7 = 1. This result rules out most of the f(R) type 
modified gravity models. However, an analysis of a La- 
grangian of the form given by Eq. (|49]) . denoted in pj} 
as f(R) = (R/a) 1+s has also been considered, and the 
conclusion is that "...this analysis is incapable of deter- 
mining whether f(R) = R 1+s gravity with 8^1 conflicts 
with Solar System tests" jl7| . Therefore, the dark matter 
model obtained from the action given by Eq. (|49| is also 
consistent with the Solar System tests of general relativ- 
ity. On the other hand, the results of also indicate 



that our model is in agreement with all observations on 
scales smaller than 10 — 20 kpc. 

Another problem facing the f(R) gravity models is the 
problem of the stability (l7L YL&i . |25| . On a time scale of 
r w 10~ 26 s a "fatal instability" develops when f"(R) < 
[18|. For the Lagrangian given by Eq. (|4"9")) we have 

f"(R) = f (v 2 g /2)(v 2 g /2 + l)^/ 2 - 1 > 0. Therefore 
this type of instability does not develop in the present 
model. 

The mass discrepancy in clusters of galaxies is a sec- 
ond major observational evidence leading to the necessity 
of considering the existence of dark matter at a galactic 
and extra-galactic scale. The total mass of a cluster of 
galaxies can be estimated in two ways. First, by tak- 
ing into account the motions of its member galaxies, the 
virial theorem provides an estimate, My. Second, the 
total baryonic mass M may be estimated by consider- 
ing the total sum of each individual member's mass. The 
mass discrepancy arises as one generally verifies that My 
is considerably greater than M, with typical values of 
My /M - 20- 30 Q. 

The virial theorem in f(R) modified gravity was de- 
rived, by using the collisionless Boltzmann equation, in 
[3ll ] . The supplementary geometric terms in the modified 
Einstein equation provide an effective contribution to the 
gravitational energy, and the total virial mass, propor- 
tional to the effective mass associated with the new geo- 
metrical term, may account for the virial theorem mass 
discrepancy in clusters of galaxies. The Lagrangian of the 
modified gravity model, which can be obtained in terms 
of quantities directly related to the physical properties 
of the clusters, and which can be determined from astro- 
physical observations, is again of the form of the action in 
Eq. (|49p . This shows that a modified gravity model, with 
an action of the form (|49|) . could give a consistent geo- 
metric description of the properties of the "dark matter" 
on both galactic and extra-galactic scales. Once the main 
physical parameters of matter at the galactic and extra- 
galactic scale, like the rotational velocities of the test 
particles around galaxies, or the intra-cluster gas tem- 
perature or the gas density profile, are known, the action 
of the modified gravity model can be completely obtained 
from observations, and the viability/non- viability of the 
model can be directly tested by using galactic and cluster 
of galaxy data, which may also offer an effective alterna- 
tive to the Solar System tests. 

Another possible physical test of the consistency of our 
model can be provided by the study of the deflection of 
light. The observational study of the propagation of light 
at the galactic or galaxy clusters level and, in particu- 
lar, the investigation of the deflection of photons passing 
through the regions where the rotation curves of mas- 
sive test particles are flat, represents one of the most 
powerful ways by which one could in principle constrain 
f(R) gravity as an alternative dark matter model for 
galactic/extra-galactic astrophysical systems. In the flat 
velocity curves region, the metric coefficients are given 
by Eqs. (TIT))) and {35]), respectively. The general analysis 
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of the lensing in the constant velocity region was per- 
formed in Ref. [H, and the obtained results can also be 
applied to the f(R) models. Therefore, lensing effects 
can in principle discriminate between the f(R) gravity 
and other dark matter or modified gravity models. By 
using the results of [3| it follows that the f(R) gravity 
model predicts slightly smaller gravitational lensing ef- 
fects in the constant velocity region, as compared to the 
standard dark matter models. 

On the other hand, the resulting deflection angle in the 
present f(R) gravity model is of the same order of mag- 
nitude as the lensing angle in the standard dark matter 
model, the isothermal sphere model, which is well con- 
firmed by observations. This shows that our solution is 
consistent with the existing observational data, and that 
high precision observations may discriminate between the 
different classes of models proposed to explain the motion 
of test particles around galaxies. 

In a series of papers [26j, using R n gravity models, a 
modified Newtonian potential of the form 



Tr 



(50) 



was considered, where m is the mass of the particle, r c 
a constant and the coefficient (3 depends on the 'slope' 
parameter n in the modified action. In the weak field 
slow motion approximation (3 can be expressed as 



Yin 2 -In- I- V36n 4 + 12n 3 - 83n 2 + 50n + 1 



6n 2 - 4n + 2 



(51) 



Using the modified Newtonian potential, given by 
Eq. (|50|) . it was found that the best fit to 15 low lu- 
minosity rotation curves in R n gravity is obtained for 
n = 3.5 26] (somewhat lower values, in particular, 
n = 2.2, were obtained in |27| and [28||, respectively). 
These results seem to suggest that a strong modification, 
i.e., a rather large value of n as opposed ton = 1, of stan- 
dard general relativity is required to explain the observed 
behavior of the galactic rotation curves. As one can see 
from Eq. (|50[) . the potential obtained in this approach is 
still asymptotically decreasing, but the corrected rotation 



curve, although not flat, is higher than the Newtonian 
one, thus offering the possibility of fitting the rotation 
curves without dark matter. However, as one can see 
from Eq. (TT?)) , the correction term to the Newtonian po- 
tential in the "dark matter" dominated region, where the 
rotation curves are strictly flat, must have a logarithmic 
dependence on the radial coordinate r. This correction 
term does not appear in Eq. (|50|) . where a power law 
modified Newtonian potential is assumed to describe the 
observed behavior of the galactic rotation curves. The 
difference is also related to the asymptotic behavior of 
the metric tensor components in the two models. These 
differences in the Newtonian limit in the two models re- 
sult in different values of the parameter n in the power- 
law modified action of the gravity. On the other hand 
we have to mention that in our approach we have com- 
pletely neglected the effect of the baryonic matter on the 
space-time geometry 

In conclusion, we have found that regions with exactly 
flat galactic rotations curves do not require any kind 
of dark matter. The rotation curves are a consequence 
of the additional geometrical structure provided by the 
modified gravity theories, and a very slight modification 
of the Einstein-Hilbert Lagrangian may account for the 
existence of "dark matter" . One can, in principle, rewrite 
the resulting field equations in terms of the Einstein ten- 
sor and interpret the remaining term as a geometrical 
energy-momentum tensor. The presence of these higher 
order terms seems to provide us with an elegant geomet- 
ric interpretation of the dark matter problem. 
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